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Abstract: Tantalizing hints of the Higgs boson of mass around 125 GeV have been 
reported at the LHC. We explore the MSSM parameter space in which the 125 GeV state 
is identified as the heavier of the CP even Higgs bosons, and study two scenarios where 
the two photon production rate can be significantly larger than the standard model (SM). 
In one scenario, T{H — )■ 77) is enhanced by a light stau contribution, while the WW* 
(ZZ*) rate stays around the SM rate. In the other scenario, T{H — )• bb) is suppressed 
and not only the 77 but also the WW* (ZZ*) rates should be enhanced. The rr rate 
can be significantly larger or smaller than the SM rate in both scenarios. Other common 
features of the scenarios include top quark decays into charged Higgs boson, single and 
pair production of all Higgs bosons in e^e~ collisions at y/s < 300 GeV. 
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1 Introduction 

Latest results from the Higgs boson search by the ATLAS [1] and the CMS [2] collaborations 
show an excess of events around the mass region of 125 GeV. The main search channel is 
the two photons decay mode of the Higgs boson, for which both experiments reported the 
rate higher than the standard model (SM) prediction. There are hints of the ZZ* decay 
mode with less significance, while no hints have been reported for the rf mode. We expect 
that the data from the current 8 TeV run will make clear the properties of the Higgs boson 
candidate. 

The Higgs sector in the minimal super symmetric standard model (MSSM) has five 
physical mass eigenstates, two CP even and one CP odd neutral scalar bosons, if CP is 
conserved in the Higgs sector, and one pair of charged scalar bosons [3]. The observed 77 
resonance at 125 GeV can be one of the three neutral Higgs bosons. Among them, the CP 
odd state [A) cannot give the 77 rate greater than that of the SM Higgs bosons, mainly 
because it lacks the W boson loop contribution to the 77 decay [4]. Among the two CP 
even Higgs bosons, both the light (h) and heavy (H) mass eigenstates can be 125 GeV and 
can have enhanced 77 rate. MSSM scenarios where the lighter of the CP even Higgs boson 
is identified as the 125 GeV state are discussed in refs.[5-21], and the possibility of the 125 
GeV state as the heavier of the CP even Higgs bosons is discussed in ref.[16]. The former 
scenario contains the so-called decoupling region where all the other Higgs bosons {H, A, 
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H^) are significantly heavier than the lighter CP even state h, whose properties resembles 
the SM Higgs boson. On the other hand, in the latter scenario where the heavier of the 
CP even state H has the mass 125 GeV, not only the mass of the lighter CP even state h 
but also those of the CP odd state A and the charged Higgs boson are bomided from 
above. 

In this study, we study carefully the consequences of this non-decoupling scenario 
of MSSM where the 125 GeV state is the heavier of the CP even Higgs bosons, H. In 
particular, we study two sub-scenarios where the two photon production rate can be larger 
than the SM. In one scenario, the — )■ 77 amplitude is enhanced by a light stau loop 
which interferes constructively with the main W boson loop, while the WW* (ZZ*) rate 
is around the SM prediction. In another scenario, the 77 rate is enhanced by suppressing 
the dominant partial decay width T{H — t- bb), and not only 77 but also WW* (ZZ*) 
production rate can be large. In both scenarios, the rf rate can be significantly larger or 
smaller than in the SM. Prediction for the mass spectra of the other Higgs bosons is also 
examined. 

The enhancement of the two photon production rate due to a light stau in the de- 
coupling region has been studied in refs.[12, 20]. We show in this report that the same 
mechanism works in the non-decoupling region as well. The suppression of T{H — )• bb) 
in the non-decoupling region has been studied in ref.[16]. We study not only the 77 and 
WW* (ZZ*) rates but also the rf rate in detail. 



2 Higgs sector in MSSM 

In this section, we briefly review the mass spectrum of the Higgs bosons in MSSM. In our 
scenarios where the two photon production rate of the heavier CP even state H is higher 
than that of the SM, relatively large Higgs couplings to the weak bosons are necessary, 
since the main contribution to the — >• 77 amplitude comes from the W boson loop [4]. 
Hence, H must be a SM-like Higgs boson. 

The MSSM Higgs sector consists of two SU (2)^ doublets, cpu and cpd which give masses 
to up type fermions and down type fermions, respectively [3]. When the electroweak 
symmetry is spontaneously broken, MSSM gives five physical mass eigenstates, two CP 
even scalar bosons h and H, one CP odd scalar boson A, and one pair of charged scalar 
bosons H^. The two CP even scalar bosons are mixed states of the real components of the 
two Higgs doublets. 



cos a — sm a 
sin a cos a 




(2.1) 



where we define h and H as the lighter and the heavier of the two CP even scalar bosons, 
respectively, whereas the current basis states and are defined as in 

Re{ct>l) = '-^^, Rei4) = '-^^. (2.2) 
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Upon the convention that the above vacuum expectation values are written as Vu = v sin /3 
and Vfi = V cos (3 with v (~ 245) GeV being the vacuum expectation value of the SM Higgs 
doublet, we can introduce another base 




sin/3 -cos/3\ / Hsm 
cos /3 sin (3 I \ H± 



(2.3) 



where Hsm is a state whose couplings to the weak bosons are the same as those of the SM 
Higgs boson, and H± is its orthogonal state which has no coupling to the weak bosons. 
From eqs.(2.1, 2.3), we have 



h = — sin (q — f3)HsM — cos (a — /3)H±, 
H = cos (a — f3)HsM — sin (a — /3)H±. 



(2.4a) 
(2.4b) 



The masses and the eigenstates of the CP even Higgs bosons in the MSSM are determined 
by diagonalizing the symmetric mass-squared matrix in the space of {H^,H^)^, 



ML Mi; 

ML Mj. 



(2.5) 



whose elements can be approximated as [22] 
3 



uu 



Mi 



1 



M2 

V2 In susy 




Mf 



12 



_3z^ 
967r2 



dd 



Ml 



327r- 



327r2 



ML- 



cos /3 



Ml + Ml + 



167r2 



3) + 



167r2 



3) 



A* ) 
(2.6a) 
(2.6b) 

32^2 ^...v^.t , ^2^2.,^ (2.6c) 

where only the leading terms for large tan (3 (tan /3 ^ 1) are kept, since large value of tan /3 
is necessary to have a SM-like Higgs boson as heavy as 125 GeV. Yj and are, respectively, 
the top and bottom Yukawa couplings in the MSSM. The soft SUSY breaking Af terms 
and the Higgsino mass /i are made dimensionless as At = At/Mgusy, Ab = A^/M, 
fit = fi/Msusy with 



rYt^Mt (i?_6)+-^y4-3A, 



'■susyi 



M. 



Mi +Mi 

tl t2 



susy 



(2.7) 



where M^^ and M^^ are masses of the stop mass eigenstates. Full analytic formulae of the 
mass matrix elements can be found in ref.[22]. By diagonalizing the matrix eq.(2.5), we 
obtain the masses of the two CP even Higgs bosons. 



Ml 



Ml,, cos^ a + ML sin^ a 



MLsm2a, 



M 



H 



M„„ sin a + M^^ cos a + sin 2a, 



(2.8a) 
(2.8b) 
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with 



sina = ^= (2.9a) 
(M^ - + (M2,)2 

cos a = ~ (2.9b) 

where we can choose the mixing angle a in the region — ^ < a < ^, since Mfj — > 
(cos a > 0) is always satisfied. The region of a can further be separated depending on the 
sign of sin a, 

1. sina<0 (M^^ < 0), -| < a < 0, (2.10a) 

2. sina>0 (M^^ > 0), < a < ^. (2.10b) 

Here the region 2 (sin a > 0) takes place if the loop contribution dominates over the 
negative definite tree-level contribution in eq.(2.6c), which can happen, for example when 
both JlAt > and > 6 are satisfied for tiny cos /3 (large tan/3). 

In the limit that the state H has the SM-like couplings to the weak bosons, we have 

|cos(a-/3)| ~ 1 (2.11) 

in eq.(2.4b). At large tan/3 (/3 ~ '7r/2), this condition eq.(2.11) selects two distinct regions, 
a — /3 ~ (a ~ '?r/2) or a — /3 ~ — vr (a ~ —-k/2). In both cases, we have cos a <^ 1, and 
the mass eigenstates in eq.(2.8) are approximated by 

~ MJ^ - 2M^^ cos a + ©(cos^ a), (2.12a) 
M|r ~ M^^ + 2M^^ cos a + ©(cos^ a). (2.12b) 

By neglecting small terms proportional to cos a, the condition that the SM-like state H is 
heavier than the other state approximately implies 

Ml-Ml>0, (2.13) 

or from eq.(2.6) 

M| + Ay^2ln^(.2y,2-M|) 

+ ^v%'A', ( 1 - + -^v%'f^\Al - p?) > Ml (2.14) 



87r2 ' ' V 12 7 327r 

Hence Ma is bounded from above by the loop contribution to the Higgs potential. 

It should also be noted from eqs.(2.6a, 2.12b) that, in order to make H a SM-like Higgs 
boson and as heavy as 125 GeV, large Msusy and ~ 6 are necessary, and we explore the 
MSSM parameter region which satisfies these conditions in the following sections. 
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3 Scenarios giving large two photon rate 



In our analysis, we consider the following Higgs production processes at the LHC, 

Gluon fusion gg ^ (p + X, (3.1a) 
Weak boson fusion qq — t- qq(j) + X, (3.1b) 
Bottom quark annihilation 66 — )• (j) + X, (3.1c) 

where (j) can be h, H or A. The SM Higgs production cross sections for the processes in 
eq.(3.1) are calculated by using the programs HIGLU [23], HAWK [24] and BBH@NNLO 
[25], respectively. The MSSM Higgs cross sections are obtained by scaling the corresponding 
SM Higgs cross sections with the ratio of the corresponding MSSM decay width over the 
SM one. The decay widths, couplings and mass spectra of the Higgs bosons and SUSY 
particles are calculated with an updated version of CPsuperH2.0 [26] which includes the 
stau contribution to the Higgs boson masses. Although the SM cross section of the bottom 
quark annihilation process is quite small compared to the dominant gluon fusion process, 
it can be significant in some MSSM scenarios. 

We consider the following constraints from the collider experiments. For the stau and 
stop masses, we adopt the lower mass bounds [27] 

Stau Mf > 81.9 GeV, (3.2a) 
Stop Mj > 92.6 GeV. (3.2b) 

Upper bounds on the e^e~ annihilation cross sections 

(T(e+e- Zh{-^ Zbb and Ztt)) , (3.2c) 
aye e~ — )• j4/i(— )• 6666, 66Tf and TfTf)) (3.2d) 

are taken from ref.[28], and those on the cross sections at the LHC 

a{pp^ h,A,H{^ Tf )) (3.2e) 

are taken from ref.[29]. Upper bound on the branching fraction 

B{t ^ bH+{^ bfi^r)) (3.2f) 

is taken from ref.[30]. Since all the physical Higgs bosons are relatively light in our non- 
decoupling scenario when < Mh ~ 125 GeV, all the above constraints in eqs.(3.2) are 
required to be satisfied in all the results presented below. In particular, significant portion 
of very large tan/3 regions is excluded by the h,A,H — )• rf and t — )• bH^{-^ fVr) search 
limits (eqs.(3.2e, 3.2f)). 

We define the ratio of a production rate at the LHC as 

p ^ <PV ^ H)B(H ^ AB) 

a{pp ^ H)SMB{H ^ AB)S^i ' ^ ' 
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which gives the H — )• AB production rate normahzed to the SM prediction. Although we 
calculate the Higgs boson cross section a{pp H + X) at = 7 TeV in this study, the 
ratio Rab should not change significantly even for ^/s = 8 TeV, since the dominance of 
the gluon fusion process remains to be valid. Hence, our results can also be applied for 
future results of ^/s = 8 TeV. Since the gluon fusion process dominates over the other 
production processes and the total decay width is dominated by T{H — )• bh) for the heavy 
CP even state H with mass around 125 GeV, Rab may be approximately written by using 
the partial decay widths, 

- (M^) ^ i^^^^ " ■ iM^) ■ (3.4) 
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By introducing a short hand notation 

_ T{H ^ ab) 

for a partial width normalized to the corresponding SM value, the production rate Rab of 
eq.(3.4) can be expressed as 

Rab ^ Tgg ■ (r^fe)"^ • r^s. (3.6) 
In this study, we examine i?^^, Rvv (y = W, Z) and Rrf, 

(3.7a) 
(3.7b) 
(3.7c) 

and identify two scenarios where the following two conditions are satisfied for the heavier 
CP even Higgs boson in the MSSM, 

123 <Mh < 127 GeV, (3.8a) 

1 <R.f.f < 3. (3.8b) 

Specifically, they are 

Light stau scenario : r^^ > 1 and rgg ■ (r^j)""^ ~ 1, (3.9a) 

Small T{H — )• 66) scenario : (r^j)"^ > 1 and rgg ■ r^j ~ 1. (3.9b) 

Since the two scenarios, light stau scenario and small T{H — t- 66) scenario, have distinct 
predictions for Rvv ^-nd Rrf, which can be tested in the current run of the LHC, we explore 
their consequences carefully in the extended parameter space of the MSSM. 
For definiteness, we explore the following MSSM parameter region, 

5 < tan/3 < 40, 110 < Mh± < 210, 
500 GeV <At< 5000 GeV, 500 GeV < ^ < 1500 GeV, 

300 GeV < Mq =Mjj = M^< 1500 GeV, (3.10a) 
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where Mjj± is the charged Higgs boson mass, Mj is the SUSY breaking sfermion mass 
parameter. The following parameters are set to fixed values, 

Ab = Ar = l TeV, 
M3 = 800 GeV, M2 = 200 GeV, Mi = 100 GeV, (3.10b) 

where Mi are gaugino mass parameters, since they do not aff'ect significantly the property 
of the Higgs bosons. The slepton soft mass parameters are explored in the region 

50 GeV <Mi = M^< 500 GeV (3.10c) 

for the light stau scenario, while it is set to a fixed value 

Mi = M^ = l TeV (3.10d) 

for the small T{H — )■ hh) scenario. 



3.1 Fermion and sfermion contributions to Vgg and r^^ 

In the SM, top quark loop contributes dominantly to the Hsm — ^ 99 amplitude. The 
bottom quark loop interferes destructively with the top quark loop for Mhsm ~ GeV, 
and for Mhsm ~ 125 GeV it counteracts the top quark loop contribution by roughly 10 %. 
In the MSSM, the heavier CP even Higgs boson, H, couples up and down type fermions, 
respectively, with the couplings 

sina / \/2m„\ tn ^^ \ 

cosQ / \f2md\ , . 

'"'^'^ = ^ ) ■ ^'-^'^^ 

Hence, when sina > in eq.(2.10b) and Mh > 30 GeV are satisfied, the bottom quark 
loop interferes destructively with the top quark loop as in the SM. When sina < in 
eq.(2.10a) and Mh ^ 30 GeV are satisfied, on the other hand, the bottom quark loop 
interferes constructively with the top quark loop, which can lead to rgg > 1. 

Similar discussion is applied for r^^. In the SM, the W boson loop contributes domi- 
nantly to the HsAi — )• 77 amplitude. The sub-dominant top quark loop interferes destruc- 
tively with the W boson loop, whereas the bottom quark loop interferes constructively, for 
Mhsm ^ 30 GeV. In the MSSM, the H coupling to the weak bosons normalized to the 
SM value is cos(a — (3). Hence, when H (~ 125 GeV) has the SM-like coupling to the 
weak bosons and tan/3 ^ 1, we have cos(a — /3) ~ sina and the top quark loop always 
interferes destructively with the W boson loop, whereas the bottom quark loop interferes 
constructively when sin a > as in the SM and destructively when sin a < 0. 

Sfermions in the third generation can have important contribution to the H ^ gg and 
the -ff — )■ 77 amplitudes due to their large Yukawa couplings. The mass eigenstates of the 
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Figure 1. r^g (solid line), r^^ (dashed line) and rgg ■ r^j (dotted line) as functions of the lighter 
stop mass for sina < (a) and sina > (b). Mh = 125 GeV, At/Mg = 2.6, Mq = = M^, 
tan 13 = 10, /i = 1 TeV, | sina| — sin/? and only the stops are considered among SUSY particles in 
the amplitudes. 



sfermions /i^2 (^^/^ < ^f2^ mixed states of the current eigenstates /l,_r with a mixing 
angle 9f, 



fA ^ cosOf -sinOA fA 
\h) ysinOf cos Of j [fRj ■ 

The mass matrix of the sfermions in the current basis is given by [31] 

/ Ml^ Ml A ^fMj^+mj + d{ mjiAf - i^rj) ^ 



(3.12) 



(3.13) 



where mj is the corresponding fermion mass and r^^ = l/r„ = tan/3 for down and up type 
fermions. The D terms are given in terms of the electric charge e/, the weak isospin /| 
and the weak mixing angle 6^ by 

dI = {if - e/ sin^ 9y,)Ml cos 2/3, (3.14a) 
D^j^ = ef sin^ O^M^ cos 2/3. (3.14b) 

The mass eigenvalues are 

Ml = "^'^ ^ "^^^ ± y{Ml, - Ml^Y + A{Ml^Y, (3.15) 
where /_ = /i and /+ = /2 are the lighter and heavier mass eigenstates, respectively, and 
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the mixing angle 9f {\0f\ < vr/2) is given by 

2mf{Af - iirf) 



sin 20 



/ 



cos 26 



h h 
M? +d{^-MI -D{ 



f 



fn 



fL 



/2 /l 



(3.16) 



(3.17) 



The heavier CP even Higgs boson, H, couples to up and down type sfermions in the mass 
eigenstate basis as follows [31] 



9h 



U±U± 



2 smQ ^ ^2 ^Qg/ _^ ON /j3 pQg2 Qr - efCos29f sin^ 9^,) 
■' sm /3 ^ V ' \ J J J J J 



^Hd±d± 



, cos a 



rrij 



=F 



V sin /3 



cos a/Li — sin aAj- sin20j, (3.18a) 



ml + M| cos(q + /3) (/| cos^ 6*/ - e/ cos 26*/ sin^ 6*^, 



cos /3 



=F 



V cos /? 



sina^u — cos aAj 



sin 261 



(3.18b) 



When is a SM-like Higgs boson with tan /3 ^ 1 and when the mixing between fi and 
fa is large, these couplings are approximated by 



9h 



^Hd±d± 



± 



2mlAl 



2m'^fj? tan /3 



■ sma. 



vcos/3(M? 



sm a. 



(3.19a) 
(3.19b) 



which are proportional to sin a, and the lighter of the mass eigenstates of the sfermions 
always interferes destructively with the top quark loop, while the heavier interferes con- 
structively with the top quark loop, independently of the sign of sin a. The lighter one 
generally contributes dominantly, and hence the squarks with large mixing always reduce 
Vgg and increase r^^ at the same time. Figure 1 shows rgg (solid line), r^^ (dashed line) 
and Vgg ■ r^^ (dotted line) as functions of the lighter stop mass for sin a < (left) and 
sina > (right), where Mh = 125 GeV, Mq = = M^, At = 2.6Mq, tan/3 = 10, 
/i = 1 TeV, I sina| = sin/3 and only the stops are considered among SUSY particles in the 
amplitudes. In both cases, the reduction of rgg due to stop contribution is always larger 
than the corresponding enhancement in r^y, and hence the light stop reduces the product 
rgg • r-y-y. When tana = — tan/3 = —10, the bottom quark contributes constructively to 
the top quark loop, giving rgg > 1 in Figure 1 (a) for large stop masses. In contrast to 
squarks, stau can increase r^^ without affecting rgg. 



3.2 Light stau scenario 

In this section, we examine R^^, Rw and Rrf in the light stau scenario of eq.(3.9a) where 
the mass of the heavier CP even Higgs boson is 125 it 2 GeV and a large R^-y is obtained 
by increasing r^^ . 
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3.2.1 R 



'77 



As discussed in Section 3.1, light stau with large mixing between tl and fji can increase 
r^-y without decreasing rgg. Hence, with light stau and heavy squarks, we can expect 



In Figure 2 (left), we show the maximum and minimum values of R-yy as functions of the 
mass difference between the lighter and heavier staus, for three different masses of the 
lighter stau, 82 GeV (sohd line), 100 GeV (dotted line) and 140 GeV (dashed hue). We 
impose in the light stau scenario the condition 



so that the enhancement of R^^ is mostly due to the light stau contribution to r^^. The 
plot shows that light stau mass difference grows, as expected, since 

large mass difference corresponds to large mixing between tl and tr. The maximum value 
of R^^ is obtained when the lighter stop has large mass ~ 1300 GeV and the stau mixing 
is large, while the minimum value of R^^ is obtained when the lighter stop has small mass 
~ 300 GeV and the stau mixing is small, within our explored parameter region of eq.(3.10), 
since light squarks generally decrease Vgg • ^-yy , as discussed in Section 3.1. 

3.2.2 Rvv 

Next we consider Ryv (^ = in the light stau scenario. Even with the help of the 

light stau contribution, significant enhancement of R^y-y over unity is possible only when the 
heavier CP even Higgs boson H has the SM-like coupling to the weak bosons as explained 
in Section 2. We therefore expect Rvv ~ ^vv ~ 1 hi this scenario, although rw < 1 
always holds. In our numerical calculation, rvv lies between 0.95 and 1, whereas Rvv is 
found to lie between 0.7 and 1.15 for the lighter stop mass between 300 and 1300 GeV, 
with little dependence on the stau masses. Rvv > 1-1 is found when sin a < where the 
bottom quark contributes constructively with the top quark loop, giving Vgg > 1.1, see 
Figure 1 (a). 

3.2.3 Rtt 

Here we discuss Rrf in the light stau scenario. As discussed above, the enhancement of 
Rjj by the light staus can be obtained for large values of fj, and tan /3 which appears in the 
Hfifi {i = —, + or 1, 2) couplings in eq.(3.19b). When /itan/3 is large, the radiative SUSY 
corrections in the bottom quark and the tau lepton masses can be important, 



^55 ■ ^^77 > 



(3.20) 




(3.21) 



Mb = -^ucos^(l + Afo) 
v2 
Y 

Mr = -^V COS /3(1 + A^) 




(3.22b) 



(3.22a) 
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Figure 2. The maximum and minimum values of R^-y (left) and Rtt (right) are plotted against the 
mass difference between the lighter and heavier staus, for three different masses of the lighter stau, 
82 GeV (solid line), 100 GeV (dotted line) and 140 GeV (dashed line) in the light stau scenario. It 
is shown that light stau mass difference grows. The maximum value of i?^-y is 

obtained when the lighter stop has large mass ~ 1300 GeV and the stau mixing is large, while the 
minimum value of R^-y is obtained when the lighter stop has small mass ~ 300 GeV and the stau 
mixing is small, within our explored parameter region of eq.(3.10), since light squarks generally 



decrease ra 



as discussed in Section 3.1. The enhanced Rtt appears when sin a < 0, while 



Rrf ~ 1 is obtained when sin a > 0, see Section 3.2.3. 



where [32] 

A{j = fj, tan /3 
A,- = tan /3 



2a 



-^MsI{M~,^ , M^^ , Ms) + Y^AJiMi^ , M,^ , ^^) 

2 2 

^MiIiMf, , Mf, , Ml) + -g^M2l(M^,, M2, /z) 



(3.23a) 
(3.23b) 



The function I(a, b, c) is given by [32] 



Iia,b,c) 



In(a7b2) + b^c^ ln(6Vc2) + c^g^ ln(cVa^ 

(a2_62)(52_c2)(a2_c2) 



(3.24) 



which is positive for all real a, b, c. The effective Higgs couplings to bb and rr normalized 
to the SM values are now given by [33] 



dHrf ■ 



cos a 
cos /? 
cos a 
cos (5 



1 



Aft 



1 + Afc 

A. 
1 + A^ 



tan a 
tan/3 
tan a 
tan/3 



(3.25a) 
(3.25b) 
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and their squared values should approximately correspond to r^g and r^f, respectively. 
The behavior of these couplings strongly depends on the mixing angle a, — ^ < a < 
(sin a < 0) in eq.(2.10a) or < a < ^ (sin a > 0) in eq.(2.10b). Figure 2 (right) shows the 
maximum and minimum values of Rrf against the mass difference between the lighter and 
heavier staus, for three different masses of the lighter stau, 82 GeV (solid line), 100 GeV 
(dotted line) and 140 GeV (dashed line) in the light stau scenario. The enhanced Rrf in 
the plot appears when sin a < and it lies roughly between 2 and 4 even when ^ 1, 
whereas Rrf ~ 1 appears when sin a > 0. Below we discuss Rrf in detail for each case. 

Rrf when sin a < 

Since H is now the SM-like Higgs boson with cos(a — /3) ~ — 1 in eq.(2.4b), we can estimate 
its deviation from the SM limit, with a small parameter e, as 

1-6. (3.26) 



tan/3 

The effective Higgs coupling to bb in eq. (3.25a) can then be expressed as 

1 



9Hbb 



1-^(2 
1 + Afe^ 



1 - e 

while the coupling to rf in eq. (3.25b) may be approximated as 



(3.27) 



9Hrf = (3.28) 

since A,- is significantly smaller than due to the electroweak couplings in eq. (3.23b). 
With the above approximation, the partial width ratio (?'bb)~^ -Trf in eq.(3.7c) is calculated 
and we find 

Rrf^r,,.( ^ \+!^' X - (3.29) 



1- Ab(l-e) 

Rrf should hence be always larger than unity if A;, is positive and rgg ~ 1. We note that 
Ab is positive in the MSSM parameter region of eq.(3.10) we explore in this study, where 
fi, M3 and At are all positive. 

Our assumption of r^g ~ 1 in the light stau scenario of eq.(3.21) can be satisfied when 
e « 2Ab/{l + 2Ab), which leads to 

Rrf ~ rgg ■ (1 + 2Abf . (3.30) 

The increase of Rrf for large mixing between fi and tr found for sin a < in Figure 
2 (right) can be explained by this mechanism, since the large mass splitting Mfj — Mf^ 
implies large ^tan/3 in eq.(3.13), which leads to large Af, from eq.(3.23a). In contrast to 
R^^, large values of Rrf are found for light squarks, because increase of A^ induced by 
light squarks in eq. (3.23a) is larger than the decrease of rgg. The maximum value of Rrf is 
obtained when mass of the lighter stop is roughly between 400 and 600 GeV and the stau 
mixing is large, while the minimum value of Rrf is obtained when mass of the lighter stop 
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is 1300 GeV and the stau mixing is small, within our explored parameter region of eq.(3.10). 



Rrf when sin a > 
In this case, we can express tana/ tan/3 (> 0) as 

tan a 



tan/3 

The effective Higgs coupling to bb in eq. (3.25a) now becomes 

1 



1-e. (3.31) 



9Hbb 



1 



1 - e 

while the coupling to rr in eq. (3.25b) may be approximated again as 



(3.32) 



1 - e 

We then find 



1 

gHrf = z • (3.33) 



2 



Our assumption of r^j ~ 1 in the light stau scenario of eq.(3.21) can be satisfied only 
when |e| ^0.1, which implies 

Rrf ^ rgg. (3.35) 

This behavior of Rrf is shown in Figure 2 (right). The plot shows little dependence of 
Rrf on the stau masses and mixing as it is expected from eq.(3.35) when sin a > 0. The 
maximum values of Rrf in the plot are obtained when mass of the lighter stop is around 
1300 GeV, while the minimum values of Rrf are obtained when mass of the lighter stop 
is between 300 and 400 GeV within our explored parameter region, since light squarks 
decrease rgg, as discussed in Section 3.1. 

Summing up this sub-section, the sub-scenarios with sin a (~ sin/3) > and sin a (~ 
— sin/3) < in eqs.(2.10) can be distinguished in the light stau scenario by measuring 
Rrf - If no significant enhancement over the SM rate is found, only the sin a > region is 
allowed, where the loop contribution reverses the sign of the off-diagonal element M^^ of 
the Higgs mass squared matrix in eq.(2.6c). 

3.3 Small T{H bb) scenario 

If the stau masses are as large as 1 TeV, their contribution to r^^ is suppressed, and the 
only alternative way to enhance R-y-y is to suppress T{H — t- bb) in eq.(3.7a), 

irbir' > 1- (3.36) 

In this section, we examine R^y-y, Rvv and Rrf in the small T{H — ?• bb) scenario of eq.(3.9b) 
where the mass of the heavier CP even Higgs boson is 125 it 2 GeV and an enhanced value 
of R^^ is obtained by decreasing r^^; . 
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Figure 3. Correlation between R^^^ and Ryv witli tlie maximum (solid line) and the minimum 
(dashed line) values of Rvv in the small T{H — ?> hh) scenario. At some point of R^^, the maximum 
value of Rvv is obtained for heavy squarks, while the minimum value is obtained for light squarks. 



3.3.1 R^^ and Rw 

When r^g is suppressed, not only R^^ but also Rvv should be enhanced. Figure 3 shows 
the correlation between R^^ and Rvv^ with the maximum (solid line) and the minimum 
(dashed line) values of Rvv for a given R^^ in the small T{H — 66) scenario. When 
we compare the maximum and minimum values of Rvv for & given R^^, the maximum 
value is obtained for heavy squarks, while the minimum value is obtained for light squarks. 
This is because light stop contributions enhance r^^ whereas their contributions to Vgg are 
common in R.^^ and Rw We find that both i?-y-y and Rvv can be as large as 3 in the 
small T[H — )• 66) scenario within our explored parameter region of eq.(3.10). 

3.3.2 Rrf 

As in the light stau scenario, Rrr in the small T{H — )• 66) scenario depends strongly on the 
sign of sin a. We show in Figure 4 the maximum and minimum values of Rrf as functions 
of the product of ^ and tan/3, for three different masses of the lighter stop, 400 GeV (solid 
line), 600 GeV (dotted line) and 1000 GeV (dashed line), when i?.y^ lies between 1.9 and 
2.1. 

The enhanced Rrf in Figure 4 appears when sin a < 0. The increase of Rrf with a rise 
in ^tan/3 is explained from eqs.(3.29, 3.23a), since Rrf increases with A;, in eq.(3.29) and 
Afe increases with /itan/3 in eq. (3.23a). Larger values of Rrf can be obtained for lighter 
stops because the reduction in r^j due to the contribution is bigger than the reduction 
in rgg, as discussed for the light stau scenario below eq.(3.30). 

The highly suppressed Rrf in Figure 4 appears when sina > 0. From eq.(3.32), r^i < 1 
can be obtained when e < 0, or tana > tan/3. The plot shows that Rrf decreases with 
a rise in /itan/3 in contrast to the case of sina < 0. This is because Rrf in eq.(3.34) 
decreases as A^ grows when e < 0. The suppression of Rrf for light stops can be explained 
as follows. First, light stops lead to small rgg, as discussed in Section 3.1. Secondly, since 
light squarks give small Vgg, r^j must get smaller in order to keep R^^ between 1.9 and 
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Figure 4. The maximum and minimum values of -R^t as functions of the product of /i and tan /3, 
for three different masses of the lighter stop, 400 GeV (solid line), 600 GeV (dotted line) and 1000 
GeV (dashed line) in the small Y(Jl — > W) scenario, when R^-^ lies between 1.9 and 2.1. The 
enhanced i?T-r appears when sin a < 0, while suppressed Rrt < 0.5 is obtained when sin a > 0, see 
Section 3.3.2. 

2.1. This requires even smaller e (< 0) in eq.(3.32), which further suppresses R^f as in 
eq.(3.34). 



4 Constraints on the other Higgs bosons 

In contrast to the decoupling scenarios of the MSSM, in our two scenarios where the heavier 
of the CP even Higgs bosons has mass around 125 GeV, and at the same time, has the 
SM-like (nearly maximum) coupling to the weak bosons, none of the other Higgs bosons 
can be very heavy. The masses of the CP odd and charged Higgs bosons are bounded from 
above when is a SM-like Higgs boson, as explained in Section 2; see eqs.(2.13, 2.14). 
From eq.(2.14), the largest Ma may be obtained when 



A, 



(4.1) 



6 is necessary from eqs.(2.6a, 



whereas in order to obtain Mh as large as 125 GeV, At 
2.12b). 

Figure 5 shows the allowed mass regions in the Mh and Ma space (lower regions) 
and in the Mij± and Ma space (upper regions) when 1.5 < -R-y-y < 2.5, for the light stau 
scenario (dashed line) and the small r(iJ — t- 66) scenario (solid line). The mass region 
with large Ma and are obtained by the large radiative SUSY correction to the Higgs 
potential in eq.(2.14). The lower bound on M}^ comes from the upper bound on the cross 
section a{e^e~ — t- Z/i), and the lower bound on M/j + Ma comes from the upper bound 
on the cross section a{e^e~ — t- AK). The main reason for the smallness of the allowed 
regions in the light stau scenario is because is constrained to be between 0.9 and 1.1 
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160 r 




Figure 5. Allowed mass regions in tlie Mh and Ma space (lower plots) and in the A/^f± and AIa 
space (upper plots) when 1.5 < R^^ < 2.5, for the light stau scenario (dashed line) and the small 
r{H bb) scenario (solid line). 

in the light stau scenario, while it is not constrained in the small T{H — ?• bb) scenario. 
Since h is non SM-like Higgs boson in our scenarios, the cross sections a{e~^e~ — t- Zh) is 
highly suppressed, while the cross sections a{e^e~ — )■ ZH) and a{e~^e~ — )■ Ah) are not 
suppressed. Hence all these Higgs bosons should be discovered in the future e^e~ collider. 
The charged Higgs boson with Mjj± < 150 GeV may be discovered from the top quark 
decay at the LHC. 

5 Conclusion 

In this study, we identify the 125 GeV state as the heavier of CP even Higgs bosons in the 
MSSM, and study two scenarios where R^^, the two photon production rate normalized to 
the SM prediction, can be significantly larger than unity. 

In one scenario, the i/ — )• 77 amplitude is enhanced by the light stau contribution 
which interferes constructively with the main W boson loop contribution. Within our 
explored parameter region, we find that R^^ as large as 2.0 can be obtained with a light 
stau near the current mass bound (Mf^ = 82 GeV), when the mixing between and tr 
proportional to /itan/3 is large and when the squarks are heavy. The WW*{ZZ*) rate, 
RvVj has little dependence on the stau masses and mixing, and we find that Rvv can be 
between 0.7 for the lighter stop mass ~ 300 GeV and 1.15 for the lighter stop mass ~ 1300 
GeV. Due to the large radiative SUSY correction to the bottom quark mass, large rf rate, 
Rtt, between 2.0 and 4.0 can be obtained even when R^^i is around unity. The maximum 
value of Rrf is obtained for the large stau mixing and for the lighter stop mass between 400 
and 600 GeV, with little dependence on the lighter stau mass. This enhanced Rrf appears 
when the mixing angle a of the CP even Higgs bosons in the basis of the two Higgs doublets 
is — 7r/2 < a < 0, while Rrf around unity is also possible when < a < -k/2. 

In another scenario, R^y-y is enhanced by suppressing the dominant partial decay width 
r{H — )• 66), and not only R^^ but also Rvv is enhanced. We find that both R^^ and Rvv 
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can be as large as 3 within our explored MSSM parameter region. As in the light stau 
scenario, Rrf can be enhanced when — 7r/2 < a < 0, while suppressed when < a < 7r/2. 
We find that Rrf can be as large as 2.5, for example, for the lighter stop mass ~ 400 GeV 
and /U tan /3 ~ 25 TeV or for the lighter stop mass ~ 600 GeV and /U tan /3 ~ 30 TeV when 
— 7r/2 < a < 0, whereas Rrf can be as small as 0.1 with little dependence on the lighter 
stop mass and on ^tan/3 when < a < tt/2, even when we assume that R^^ lies between 
1.9 and 2.1. 

We also study mass spectra of other three Higgs bosons, h, A and in both scenarios 
when R^^ lies between 1.5 and 2.5. We find in both scenarios that all the masses are 
bounded from above within our explored parameter regions, such that single and pair 
production of all the Higgs bosons should be observed in e^e~ collisions at ^/s < 300 GeV. 
The charged Higgs boson mass should lie in the region 120 < Mjj± < 150 GeV and may 
be discovered in the top quark decays at the LHC. 
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